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ABSTRACT 

A standard method to study the mass distribution in galaxy chisters is through 
strong lensing of background galaxies in which the positions of multiple images of the 
same source constrain the surface mass distribution of the cluster. However, current 
parametrized mass models can often only reproduce the observed positions to within 
one or a few arcsec which is worse than the positional measurement uncertainty. One 
suggested explanation for this discrepancy is the additional perturbations of the path 
of the light ray caused by matter density fluctuations along the line of sight. We inves- 
tigate this by calculating the statistical expectation value for the angular deflections 
caused by density fluctuations, which can be done given the matter power spectrum. 
We find that density fluctuations can, indeed, produce deflections of a few arcsec. We 
also find that the deflection angle of a particular image is expected to increase with 
source redshift and with the angular distance on the sky to the lens. Since the light 
rays of neighbouring images pass through much the same density fluctuations, it turns 
out that the images' expected deflection angles can be highly correlated. This implies 
that line-of-sight density fluctuations are a significant and possibly dominant system- 
atic for strong lensing mass modeling and set a lower limit to how well a cluster mass 
model can be expected to replicate the observed image positions. We discuss how the 
deflections and correlations should explicitly be taken into account in the mass model 
fitting procedure. 

Key words: gravitational lensing: strong - gravitational lensing: weak - galaxies: 
clusters: general 



> 1 INTRODUCTION 



X 



Galaxy clusters remain important tools to study the uni- 
verse as a whole as well as interesting objects in their own 
right. Their rarity and recent formation means they are es- 
sential for understanding the formation of structure in the 
Universe and their high mass implies that they are bright 
X-ray sources and act as gravitational lenses. Both X-ray 
observations and weak and strong lensing can be used to in- 
fer the properties of the cluster mass distribution. A^-body 
simulations of structure formation tell us that relaxed clus- 
ters should have a universal mass distribution which can be 



described succinctly by the NFW profile (Navarro, Frenk 



& White 19971 or similar parametrizations. Since clusters 



formed late, they are predicted to be less centrally concen- 



trated than the mass profile of galaxies or groups (Bullock 
et al.|2001| ) 



Strong lensing provides information about the surface 
mass distribution in the innermost region of a galaxy clus- 



ter. That region is critical for measuring the inner slope and 
the concentration parameter of the mass profile, which can 
be compared to the predictions of simulations. For exam- 
ple, we can take a specified model for the lensing cluster 
mass distribution and fit the parameters of that model so as 
to reproduce the observed positions of the images with the 
least scatter. This is a common approach in the literature 
and there are publicly available codes such as GRAVLENS 



(Keeton 20011 and lenstool (Julio et al. 20071 for that 



purpose. The NFW profile has been found to be a reason- 
able model for several galaxy clusters although there are also 
cases where a shallower inner slope is preferred. However, in 
many cases the best-fit mass model can only reproduce the 
observed image positions to within a few arcsec, particu- 
larly when there are several sources at a range of redshifts. 
Such precision is unsatisfactory given the sub-arcsec preci- 
sion of the image positions in HST imaging. For example, 
for the much studied cluster Abell 1689, |Broadhurst et aL] 
( |2005[) found a root-mean-squa re offset per image of 3.2 arc- 
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Halkola, Seitz & Pannella (20061 found 2.7 arcsec, and 



Limousin et al. (20071 reported 2.9 arcsec. Other examples 
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are Abell 38 3 (1.95 arcsec, [Zitrin et al.|[20TT| , Abell 2218 
(1.4 9 arcsec, |Eli'asd6ttir et al.||2007[ ), Abell 1703 (1.4 arc- 
sec,|Limousin et al.|2008| ), and CL0024+ 1654 (~ 2.5 arcsec, 
Zitrin et al.|2009| ) 



A possible explanation for these offsets is that the path 
of the image light ray is not only strongly lensed by the 
galaxy cluster but also weakly lensed by matter density fluc- 
tuations all along the path. This means the size of the offsets 
will depend on the distribution and growth of structure be- 
tween the source redshift and now. We refer to this process 
as cosmic weak lensing (CWL), to avoid confusion with the 
weak lensing signal from the galaxy cluster itself. We inves- 
tigate the impact of CWL by calculating the expected rela- 
tive image deflection angle produced by matter fluctuations 
from the matter power spectrum, using a formalism from 
lensing of the cosmic microwave background. We examine 
the dependence of the relative deflection on source redshift 
and on the image position relative to the cluster and the 
source. Finally we discuss how the line-of-sight matter fluc- 
tuations can be taken into account when fltting parametric 
mass models. 

We work in the Born approximation, evaluating the 
lensing effects along the unperturbed line of sight and we 
assume flat space for simplicity. We make use of the re- 
sults and notation of the CMB lensing review by [Lewis fc] 
Challinor ( 2006 1 . For the numerical examples we assume a 
standard flat ACDM cosmology with h = 0.71, Qb = 0.045, 
Qm = 0.27 and Qa ~ 0.73, which is very similar to the 



best-fit model of the WMAP 7- year data (Komatsu et al 



20111. We use the iCoSMO package (Refregier et al. 20111 
to calculate the evolution of cosmological quantities. 



2 METHODOLOGY 

We want to calculate the deflection angle due to CWL of an 
image which is strongly lensed by a galaxy cluster. Concep- 
tually, the observed image position 6 is the sum of the source 
position /3, the strong lensing deflection caused by the clus- 
ter asL and the weak lensing caused by cosmic weak lensing 
a CWL, 



f3 ^ - oisL — acwL 



(1) 



We cannot expect to calculate acwL for a speciflc image. 
Instead, we can calculate its expectation value by treating 
the path of the image light ray as a random line of sight. 
This means that the statistics of the deflection angle are de- 
termined by the statistics of the density fluctuations, i.e. by 
the 3D matter power spectrum. 

The CWL deflection angle of image i is related to the 
lensing potential sourced by the matter along the line of 
sight through oci — Vi/j^. The angular power spectrum of the 
lensing potential is {tpi,%l)j), where the potential is projected 
from the sources at comoving distances Xi s^nd Xj- In the 
Limber approximation {£ 3> 1) the lensing potential power 
spectrum can be written as 



8^ 

f3 



xdxr^Wx,x) [ 

V XiX 



Xj - X 
XjX 



(2) 



As usual in the Limber approximation, only the transverse 
modes contribute to the power spectrum so the integration 
extends up to the closest source, x* = niin(xi, Xj)- "P^ik, x) 



is the 3D power spectrum of the Weyl potential. This can 
be related to the matter power spectrum P{k, x) through 



V^ik,x) = 



9nl{x)H\x) P{k,x) 



(3) 



87r2 k ' 

where Qrn is the matter density and H is the Hubble rate. 

Now we can calculate the variance of the deflection an- 
gle fleld from the angular power spectrum of the potential 
(see Lewis fc Challinor|2006 for details), 

C,i{r) = (4) 
= (VV.-VV-,) (5) 

= ^^JmJ xdxr^Wx,x)F^j{x)Jomx)), (6) 

where r{x) is the angular separation between the light rays 
producing the images (which will vary with distance for 
strongly lensed images) and Fij (x) is shorthand for the lens- 
ing efficiency factors in equation ([2|. 

There can be no observable consequences if CWL causes 
the whole observed field to be displaced coherently. Hence, 
the relevant issue to consider is not how large the deflection 
angle Ui of image i is, but rather how large is the difference 
of cti and the deflection a' of a fiducial line of sight. This 
reference sight-line determines a 'zero-point' for the CWL 
defiection angles in the observed fleld. It can be arbitrarily 
chosen in principle, but it is convenient to choose one that 
passes through the centre of the cluster lens model. The 
variance of the relative deflection is given by 

al = {{a, - a'f) ^ 2CgiiO) - 2Cg,ir), (7) 

where for both terms the integral over x in equation ([6| 
should be evaluated up to the source distance Xi- This is the 
expected offset of an image relative to the cluster centre. The 
strong lensing event causes the angular separation to vary 
with comoving distance. 



■ rix) 



ro for < X < Xd, 
(ro - rs)^f5^ for Xd < X < Xi 



and rs is the separation in the image plane between the 
source and the reference position. The covariance of two im- 
ages of different, or identical, sources is 



{{a^ - a') ■ (aj - a')) 



(8) 



where again the angular separations will vary with comoving 
distance and all terms are evaluated at x* = niin(xi,Xi)- 

How does this relate to the cluster mass modeling? Let 
us summarize how parametric modeling proceeds: the best 
model is the one that reproduces the observed positions {/3j} 
with the least amount of scatter and is typically found by 
minimizing a x^-function 

- IJ'i,s)^Tj^Wj,s - fJ.j,s), (9) 



x' = ^(ft,. 



where 5 labels the RA and Dec components, the model- 
predicted position of image i is /i^ and C is the image posi- 
tion covariance. The number of degrees of freedom for Mm 
image constraints of A^so sources and a model containing Np 
parameters is 2{Nim — Nso) — Np since each observed position 
consists of two coordinates and the source position of each 
image family must be inferred as well. 
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Figure 1. Expectation of the relative deflection an of an image 
as a function of distance on the sky to the cluster centre, for 
various source redshifts Zs- In this case the strong lensing event 
has been disregarded so that r(x) = t. Hence, this figure shows 
the typical scale of the relative deflection but it is not exact for 
a strongly lensed image. 



In the literature, the uncertainty on the image position 
is most often assumed to be the positional measurement un- 
certainty which is of the form cr^j = a Sij, i.e. the covariance 
matrix is diagonal and the diagonal elements are the same. 
However, the CWL covariance matrix is more complex with 
diagonal elements af^ and off-diagonal elements crfj , as given 
by equations ^ and ([8|. Hence, using the CWL covariance 
matrix will formally lead to a best-fit model which is dif- 
ferent from the one found when only taking the positional 
measurement uncertainty into account. 

We have stipulated that we treat the image sight-lines 
as random, in the sense that no other information than the 
source redshift and position is being used. In practice, this 
means that we can use the standard 3D matter power spec- 
trum P{k, x)- A more refined treatment would also take into 
account the presence of the lensing cluster, which is a rare 
type of object located at the intersection of filaments in the 
cosmic structure. Hence there are likely more extreme den- 
sity fluctuations in the surroundings of the cluster than what 
we assume. This could perhaps be modeled by biasing the 
input matter power spectrum in a suitably calibrated fash- 
ion. We leave this to future work. 

In short, the expected deflection angle of an image rela- 
tive to the cluster centre is given by an, equation ([7|. Given 
a set of observed multiple image positions, the variances 
together with the covariances afj (equation (j8|) form the 
CWL covariance matrix which should be used to find the 
best-fit mass model. 



3 NUMERICAL EXAMPLES 

Let us now discuss two examples to illustrate typical results. 

First we consider a single image located some distance 
from the cluster centre, and for simplicity we disregard the 
strong lensing event, i.e. the angular separation between the 
unperturbed light ray and the sight-line through the clus- 



ter is constant. Fig. [T] shows the resulting expected relative 
defiection an and how it increases with the separation and 
source redshift. It is clear that CWL defiections can easily 
produce defiections of a few arcsec, increasing with separa- 
tion on the sky, and also that a large fraction of the total 
defiection is contributed from evolved structure in the late 
universe. Fig. [l] can be used to estimate the scale of the rel- 
ative deflection for a given redshift, but for strongly lensed 
images the specific r(x) of the path of the light ray must be 
taken into account for a precise result. 



The second example is Abell 1689 (e.g. Broadhurst et al 



(20051; Limousin et al. (2007 1 ) which has a very large num- 



ber of multiple images. We consider a small subset of these 
to illustrate the calculations. Source 2 at 2 = 2.5 has five 
images located between 6 and 60 arcsec from the centre of 
the cluster while source 4 at 2 = 1.1 also has five images, 
of which we consider two that are both ~ 40 arcsec from 
the centre. Fig. [2] shows the structure of the held. We cal- 
culate all the elements of the covariance matrix a^j, taking 
the varying separations r{x) into account by assuming rea- 
sonable choices for the true source positions. The expected 
relative deflection angles an of these images w.r.t. the cluster 
centre are given in Table [l] and shown in Fig.|2]as error cir- 
cles. Again, it is clear that the impact of CWL increases with 
the angular separation of the images and with the redshift 
of the sources. E.g., images 4a and 4b have expected rela- 
tive deflections which are half as large as that of image 2c, 
even though all three images are about equally far from the 
cluster. Another case is image 2e which has a much smaller 
expected deflection than 2a-2e since it is very close to the 
cluster centre. In Limousin et al. ( 2007 1, the rms image plane 



offset per image for system 2 was 1.62 arcsec. Taking a naive 
rms average of the flve an (and hence disregarding correla- 
tions) we obtain an expected offset per image of 1.52 arcsec, 
i.e. a very similar number. The correlation matrix (see Table 
[1]) shows various degrees of correlation between the deflec- 
tions of the images, e.g. it can be seen that images 2b and 
2c which are relatively close together have strongly corre- 
lated deflections with a correlation coefficient of 0.88. This 
is of course not surprising since the light rays producing 
each image must have passed through much the same den- 
sity fiuctuations. On the other hand, images 2d and 4b are 
even closer together but have a slightly smaller correlation 
coefficient of 0.82 due to the lower redshift of source 4. 



4 APPLICATIONS 

While CWL can explain at least a significant part of the off- 
sets between best-flt model prediction and observed image 
position, it should really be taken into account already at 
the point of model fitting. As discussed above, the covari- 
ance matrix afj can be inserted straight-forwardly into exist- 
ing methods based on x^-minimization in the image plane. 
In practical strong lensing applications, it is not necessarily 
feasible to treat the source positions as free parameters to 
be marginalized over. Instead, the family of images are pro- 
jected backwards and the mean of these images is treated 
as a proxy for the source position which can then be pro- 
jected forwards to the image plane. While this is not prob- 
abilistically rigorous, it reduces the computational load sig- 
nificantly. It should be relatively straight-forward to adapt 
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Table 1. Multiple images lensed by Abell 1689: source redshift, observed distance on the sky to the 
cluster centre, expected relative deflection an, and correlation matrix . 



ID 


z 


r 1 arcsec 


an 1 arcsec 


Corr 












2a 


2.5 


61.9 


1.97 


1.00 












2b 




48.7 


1.84 


-0.29 


1.00 










2c 




42.6 


1.65 


-0.03 


0.88 


1.00 








2d 




32.4 


1.17 


-0.45 


0.11 


-0.20 


1.00 






2e 




6.0 


0.40 


-0.48 


0.53 


0.32 


0.60 


1.00 




4a 


1.1 


43.3 


0.89 


-0.55 


0.59 


0.30 


0.53 


0.48 


1.00 


4b 




40.0 


0.82 


-0.48 


0.20 


-0.10 


0.82 


0.49 


0.78 
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Figure 2. Schematic representation of the central 2x2 arcmin 
of the Abell 1689 field showing the five images (a-e) of the back- 
ground source 2 at z = 2.5 (following the usual labeling in the 
literature) as well as two of the five images of source 4 (z = 1.1). 
The centre of the cluster is marked with a '-|-'. The red circles 
show the Ian error circles as given by equation ([tJ, assuming 
that the true position of the sources are at the marked positions. 
Table^gives the numerical values as well as the image correlation 
matrix. 



our treatment of LSS to such a framework. Each backwards 
projection could be weighted by al^ the resulting weighted 
mean would be the assumed source position, and the rest of 
the analysis could proceed as discussed above. There are also 
methods in which the x^-function is defined in the source 
plane, which is computationally even more efficient but suf- 
fers from possible biases. It appears to be more of a challenge 
to use our method in such a case, since the covariance ma- 
trix alj would need to be transformed back to the source 
plane while taking into account both shear and magnifica- 
tion effects. 

A different issue relates to the small size of strong lens- 
ing fields which is between some tens of arcsec across and 
one or two arcmin. This means we are predicting CWL ef- 
fects on extremely small scales and we must worry about 
our ability to calculate P(k,x). Fig. [s] displays the inte- 
grand of Cgii^') — Cgi{r) and it can be seen that most of 



the signal comes from the region of parameter space where 
fc = 1 — 10/i/Mpc. This is far into the nonlinear region 
where baryonic effects also play a very significant role. The 
results discussed above were calculated using the popular 



halo model prescription of Smith et al. ( 2003 1 which is cali- 



brated on dark matter-only A'^-body simulations and cannot 
be expected to deliver precise power spectra for very high 
k. There are simulations of structure formation including 
baryonic processes which are improving rapidly but there 
does not appear to be consensus regarding the various feed- 
back processes that play an important role at the moment. 
To estimate the possible impact of this issue we have used 



the recent simulation results of van Daalen et al. (20111, 
who measured the small scale matter power spectrum us- 
ing the OWLS suite of simulations ( jSchaye et al.|[2010[ ). 
Their 'most realistic' simulation based on the WMAP7 cos- 
mology follows the formation of structure down to scales 
of ~ 500/i/Mpc and it includes AGN feedback and does 
not suffer from the over-cooling problem. We take the mat- 
ter power spectrum, which has about twice as much power 
at fe = 100/i/Mpc, from that simulation for the range 
fc > 1 ft/Mpc for all redshifts as an input and test how much 
the example results change. Comparing the case of parallel 
light-rays (Fig.[T|, for a source redshift of 2s = 1 the relative 
deflection angle an increases by 9% for a separation of r = 3 
arcsec and by 2% for r = 30 arcsec, while for a source red- 
shift of Zs — 7 the corresponding increases are 11% and 1%, 
respectively. Hence, the effect of the changed matter power 
spectrum is relatively benign for large angular separations 
while for small angular scales the difference is smaller than 
a tenth of an arcsec, and hence less than the pixel size in 
HST imaging. Hence, we conclude that uncertainty about 
the small scale power spectrum is a minor complication at 
present. 



5 DISCUSSION 

Finally, we compare our analysis with the literature and 
summarise. 



D'Aloisio & Natarajan (20111 investigated the effect 



of CWL on the positions of strongly lensed images by 
ray-tracing through multiple lens planes, populated with 
NFW halos from the Millennium A^-body simulation (MS, 



Springel, Frenk fc White (20061). They found that the mean 



displacement of the light-rays was 9kpc (comoving), corre- 
sponding to 2.3 arcsec at z = 0.2. Since they average over a 
distribution of source redshifts it is difficult to make a direct 
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Figure 3. Where does the relative deflection angle signal come 
from? The contours show the integrand of o-^^ = 2Cgi (0) — 2Cgi (r) 
for r = 10 arcsec (see equation(j6]l). The lines show the co-moving 
wavenumber k of the matter power spectrum. Most of the sig- 
nal contributing to cr?^ comes from the region between 1 and 
lOh/Mpc. For smaller r, the signal region will move to the right. 



comparison with our results, but we do note that their re- 
sults are comparable in order of magnitude. They also find 
that structure associated with the cluster lens ('correlated 
large-scale structure') can perturb the image positions, but 
to a lesser degree than CWL. In a similar analysis also based 



on the MS, Julio et al. ( 2010 1 found that CWL added deflec- 



tions of typically 1 arcsec while correlated large-scale struc- 
ture contributed 0.6 arcsec. They also noted that, in a few 
cases, CWL could change the multiplicity of images, an ef- 
fect which we do not take into account here. 

Unlike for the strong lensing regime, the impact of CWL 
on the weak lensing signal measured further away from the 
centre of a cluster lens has been studied in detail, e.g.|Hoek- 



stra 



(20031 or Dodelson (20041. Recently, Becker fc Kravtsov 



(20111 and Hoekstra et al. (20111 confirmed in simulations 



that CWL is an irreducible systematic to WL cluster mass 
measurements, but that it can be accounted for stochasti- 
cally, i.e. in a similar spirit to the present work (see also 
Gruen et aT] ( |2011[ )). In [Umetsu et"aL| \2Qll\ and |Zitrin| 
et al. (20111, the CWL contribution to the measured sur- 



face mass density in radial bins was taken into account for 
a joint strong lensing-weak lensing analysis, using the pre- 
dicted bin-to-bin covariance calculated from cosmic shear 
power spectra. 

Another proposed explanation for the positional off- 
sets is substructure in the lensing cluster. Individual clus- 
ter galaxies can be modelled concisely but coarsely through 



scaling relations, and Julio et al. (20101 found that a 20% 



random variation in the scaling relations could introduce 
scatter in the image positions of up to 1 arcsec. This ap- 
pears to be slightly less significant than our prediction for 
the effect of CWL, but we do not claim that CWL is the 
dominant source of image position uncertainty, only that it 
is a likely explanation for the reported offsets. 

In summary, we have presented a reasonable explana- 
tion for the offsets between the observed and the model- 
predicted image positions which are often encountered in 



current parametric modeling of strong lensing galaxy clus- 
ters. We find that additional lensing by density fluctuations 
along the line of sight - cosmic weak lensing - can displace 
images by a few arcsec, which is similar to the discrepan- 
cies reported in the literature. Since CWL is an unavoidable 
and possibly dominant systematic, we propose to account 
for this by fitting the mass model according to the image 
position covariance matrix of CWL. We intend to do this 



for the CLASH sample (Postman et al. 20111 of 25 galaxy 



clusters, which are currently being observed. This will allow 
us to infer parametric mass models with much better under- 
standing of the uncertainties on the model parameters. 
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